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Abstract. 

Exactly-solvable model of the linear singular oscillator in the relativistic 
CO ■ configurational space is considered. We have found wavefunctions and energy spectrum 

| for the model under study. It is shown that they have correct non- relativistic limits. 
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1. Introduction 



In the framework of the non-relativistic quantum mechanics Refs. [1-5] studied in detail 
the linear singular oscillator, which is described in the configuration representation by 
the Hamiltonian 



h 2 d 2 muj 2 2 9 /, ,n 

The singular oscillator, being one of the rare exactly solvable problems in the non- 
relativistic quantum mechanics has been extensively used in many applications. There 
are a lot of quantum mechanical and quantum field theory problems leading to the 
solution of the Schrodinger equation with Hamiltonian (jl.lj) . For example, it served as 
an initial point in constructing exactly solvable models of interacting iV-body systems 
[2j. It was also used for the modeling of the diatomic [H] and polyatomic [Zj molecules. 
The interest to the given Hamiltonian recently has increased in connection with its use 
for the description of spin chains [S], quantum Hall effect jU], fractional statistics and 
anyons [10J. 
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In the case of a constant frequency oj, a complete set of orthonormalized 
eigenfunctions of Hamiltonian (jl.ljl in the interval < x < oo can be choosen in 
the form (see, for example [3]) 



C nrci (^) 



/mu\ d+l n\ 
\~h) T(d + n + l) 



i 



x d + 2e 2n x L d n \^—x z J , (1.2) 



where d = lyl + ^p 2 ) an d -^(y) are the associated Laguerre polynomials. The 
corresponding eigenvalues of if at have the form 

E nonrel = + rf + l), 71 = 0, 1, 2, 3, ... . (1.3) 

The purpose of this paper is to construct and investigate a relativistic exactly 
solvable model of the quantum linear singular oscillator (jl.ljl . which can be applied for 
study relativistic physical systems as well as systems on a lattice. 

Our construction is based on the version of the relativistic quantum mechanics, 
which was developed in several papers and applied for solution of a lot of problems in 
particle physics [11-24]. Although this version of the relativistic quantum mechanics is 
closely analogous to the non-relativistic quantum mechanics, its essential characteristic 
is that the relative motion wavefunction satisfies a finite-difference equation with the 
step equal to the Compton wavelength of the particle, A = h/mc. For example, in the 
case of a local quasipotential of interaction V (r) the equation for the wavefunction of 
two scalar particles with equal masses has the form 

[H + V (r)] i>{r) = E^{r) , (1.4) 

where the finite-difference operator Hq is a relativistic free Hamiltonian 



Hq = mc 



iX L 2 

cosh (iXd r ) H sinh (iXd r ) H k exp (iXd r ) 

r 2 (mcr) 



1.5) 



and L is the square of the angular momentum operator and d r = j^. The technique 
of difference differentiation was developed and analogues of the important functions of 
the continuous analysis were obtained to fit the relativistic quantum mechanics, based 
on Eq. ((131) H2H3]. 

Unlike non-relativistic quantum mechanics, because of presence of the finite- 
difference operator e lXdr in presented relativistic quantum mechanics both the 
Hamiltonian and the quasipotential contain the imaginary terms. However, they are 
hermitian operators with the respect to scalar product j ip*ipdr, where ip (f) and ip (r) 
are the square-integrable functions. 

We note that, there is a regular method for construction of the quasipotential in 
the framework of the field-theoretical formalism. However, it can be also introduced 
phenomenologically. 

The space of vectors r is called the relativistic configurational space or f-space. 
The concept of the r-space has been introduced for the first time in the context of the 
quasipotential approach to the relativistic two-body problem 
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The quasipotential equations for the relativistic scattering amplitude and the wave 
function ip (p) in the momentum space have the form [To"! l2"o] 

A (p, a) = |V (p, q- E q ) + -JL J V (p, k- £ ff ) G q (k) A (k, g) aM k , (1.6) 
i> (p) = (2vr) 3 S (p (-) q) + -^G q (p) J V (p, k; E q ) V (£) d«k , (1-7) 



n2 



where 



and V (j), k; E q j is the quasipotential. 

The integration in (jl.fij) and (jl.7j) is carried over the mass shell of the particle with 
mass m, i.e. over the upper sheet of the hyperboloid po 2 — p 1 = m 2 c 2 , which from the 
geometrical point of view realizes the three-dimensional Lobachevsky space. The group 
of motions of this space is the Lorentz group SO (3, 1). 

Equations (jl.fij) and (jl.7j) have the absolute character with respect to the geometry 
of the momentum space, i.e., formally they don't differ from the non-relativistic 
Lippmann-Schwinger and Schrodinger equations. We can derive Eqs. (jl.fij) and (jl.7j) 
substituting the relativistic (non- Euclidean) expressions for the energy, volume element, 
and ^-function by their non-relativistic (Euclidean) analogues: 

„2 

E q = ^--*E q = ^q*c 2 + m 2 c 4 , 
2m 

dk 

dk -> dn k = , (i.9) 

J- I m .2..L> 



5 (p - q) -> 5 (p(-)q) . 

As a consequence of this geometrical treatment, the application of the Fourier 
transformation on the Lorentz group becomes natural instead of usual one. In this case 
the relativistic configurational r-space conseption arises. 

Transition to relativistic configurational r-representation 

V> (f) = — ^ [ £ (p, f) V (p) dQ p (1.10) 

is performed by the use of expansion on the matrix elements of the principal series of 
the unitary irreducible representations of the Lorentz group: 
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l-tr/A 

Po — pn^ 



\ mc / 
r = rn, < r < oo, 

ft = (sin 6* cos cp, sin 6 1 sin cp, cos 9) , p = \J p 2 + m 2 c 2 . 

The quantity r is relativistic invariant and is connected with the eigenvalues of the 
Casimir operator C = N 2 — L 2 in the following way 

C = X 2 + r 2 , (1.12) 

— * — * 

where L and N are the rotation and boost generators. 

It is easy to verify that the function (the relativistic "plane wave") (jl.llj) obeys the 
finite-difference Schrodinger equation 

= 0. (1.13) 

The relativistic plane waves form a complete and orthogonal system of functions in 
the momentum Lobachevsky space j26| . 

If we perform the relativistic Fourier transformation (jl.K)Jl in Eq. (jl.7j) . we arrive to 
the finite-difference Schrodinger equation (jl.4j) with the local (in general case, non-local) 
potential in the relativistic f space. 

In the relativistic r space the Euclidean geometry is realized and, in particular, 
there exist a momentum operator in the relativistic configurational f representation 



p=-fi (e iXdr - H ) - m-e iXdr , (1.14) 
where a three-dimensional vector rh has the following components 

d sirup d 



mi = i yzos <p cos 9 
mi = i I sin <p cos " 



39 sin 9 dip 
d cosip d 
89 sin^ dip 



d 

m 3 = -zsm9—. (1.15) 
69 

The components of ()1.14j) and free Hamiltonian obey the following commutation 
relations: 

\PhPj] = [Pi, H ] = 0, i,j = 1, 2, 3 . 

The relativistic plane wave is the eigenfunction of the operator p: 

p£(p,r) =pXip,r) . (1.16) 

This means that (jl.llj) describes the free relativistic motion with definite energy 
and momentum. 
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In the non-relativistic limit we come to the usual three-dimensional configurational 
space and relativistic plane wave goes over into the Euclidean plane wave, i.e. 

lim £ (p, r) = e^ h . (1.17) 

Note that all the important exactly solvable cases of non-relativistic quantum 
mechanics (potential well, Coulomb potential, harmonic oscillator etc.) are exactly 
solvable for the case of Eq. (jl.4j) . too. 

2. Relativistic Quantum Mechanics: The One-Dimensional Case 

In the one-dimensional case the relativistic plane wave takes a form [TH] 

/ \ —ix/\ / \ ix/X 



mc J \ mc 

or, in hyperpolar coordinates 

Po = mccoshx, p = mcsinhx, (2.2) 

we have 

ap,x) = e l ^\ (2.3) 

where x = hi (^^r) * s rapidity. 

The one-dimensional plane waves obey the completeness and orthogonality 
conditions 



oo 

1 



2nh 

— oo 
oo 

1 



2nh 



£ (p, x ) C (Pi x ') d£l p — 5 (x — x') , (2.4a) 
f (p, x ) € (p', x)dx = 5 (p(-)p') = 5 (mc (x ~ x')) , (2-46) 



where dVt p = rnc^ = mcdx is the invariant volume element in the one- dimensional 
Lobachevsky space, realized on the upper sheet of the hyperbola p\ — p 2 = m 2 c 2 , p > 0. 
The free Hamiltonian and momentum operators are finite-difference operators 

H = mc 2 cosh (i\d x ) , p = — mcsinh (i\d x ) . (2.5) 

Plane wave (j2.1|) obeys the free relativistic finite-difference Schrodinger equation 

(h - £ p ) £ (p,x) = , E p = cp Q = c^p 2 + m?c 2 , (2.6) 

and 

Pi (P, x) = pi (p, x) . 
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3. The Finite-Difference Relativistic Model of the Linear Singular 
Oscillator 

We consider a model of the relativistic linear singular oscillator, which corresponds to 
the following interaction potential 



V{x) 



—muo 2 x (x + iX) + "' 



e iAa * , (3.1) 



2 x (x + iX) 

where g is a real quantity (we will assign a restriction for the values of the parameter g 
below). 

For g = it coincides with quasipotential of the relativistic linear oscillator studied 
in detail in [TH] . 

Let us note that in contrast to the case of the Coulomb potential EH1 121] , which 
can be calculated as an input of the one-photon exchange, the relativistic generalization 
of the oscillator or singular oscillator potential is not uniquely defined. Therefore, for 
construction of the quasipotential (j3.1|) we proceed from the following requirements for 
the quasipotential: a) exact solubility; b) the correct non-relativistic limit; c) exsistence 
of the dynamical symmetry. 

The operator (j3.1|) is hermitian with respect to a scalar product 

oo 

frM = J ip*(x)<p(x)dx , (3.2) 

— oo 

i.e. (Vip,(p) = (ij},V(p). Here functions ip{x) and (p(x) vanish at x — ±oo together 
with all their derivatives. In this connection, we note that, the hermitian conjugate of 
operator A = f(x)e ldx with the respect to the scalar product (|3.2j) has a following form: 

A+ = e -"i* dx f*{x) , 

where f(x) is some complex function. 

It is to be emphasized that potential (j3.1|) possesses the correct non-relativistic 
limit, i.e. 

lim V(x) = -mu 2 x 2 + -^r . 

c-oo V ' 2 X 2 

A relativistic singular oscillator is described by the following finite-difference 
equation 



mc 2 coshiAda + -muo 2 x(x + iX)e iX9x H — - — — e iXdx 

2 x (x + iX) 



i/){x) = Eif>(x), (3.3) 



with the boundary conditions for the wavefunction ip(0) = and ip(oc) = 0. 
We shall confine ourselves to the interval < x < oo. 

In in terms of dimensionless variable p — x and parameters ujq = go = j§- the 
equation ()3.3|) takes a form: 



cosh^ + i^e^ + ^e^ 



1>(p) = —MP) , (3-4) 



mc 2 
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where = p(p + i). 

To solve the equation (|3.4|) we choose ip(p) as 

= c(-p) WwjfT (1/ + ip) = c{-p)^M{p)Q{p) , (3.5) 

where a and z/ are arbitrary constant parameters, which will be defined below. 
Functions 

(-P) ia)=ia T Y P (t) a) M(P) = ^ ( " + ' P) (3 ' 6) 

are connected with the behaviour of the wavefunction ^>(p) a ^ points p = and p = oo, 
respectively. 

Inserting ()3.5|) into (|3.4|) . we obtain: 

fi(p) = 2iepQ{p) , (3.7) 



^o 2 (P (2) ) +P (2) +2^o ^ 



(a + ip) (u + ip) e " 5- —7 : r^-e 

where e = E/hv is dimensionless energy. 

Now we choose constant parameters a and z/ in such a way, that they satisfy the 
following relation: 

*W + -(«-»)(•>-»)■ (3-8) 



uJq (a - 1 + ip) {y — 1 + *p) 
Hence we get 



11/ 2 



« = g + 2 V 1 + ^ ( 1 " V 1 - %o^o 2 ) , (3.9a) 



11 2 



17 = 2 + 2 V 1 + ^ V 1 + V 1 " %0 ^J • (3 ' 96) 
Then the function fi(p) will satisfy the difference equation 

[(a + ip) [y + ip) e~ idp - (a - ip) (u - ip) e idp ] Q(p) = 2iepVl{p) . (3.10) 

By substitution of the f2(p) function expansion as following power series 

00 

n(p) = X>(ip)* (3.11) 

fc=0 

into ()3.1()|) one finds that the coefficients at odd degrees of ip become zero (i.e., all 
e2fc+i = 0, k = 0, 1, 2, 3, . . .), but the coefficients at even degrees of ip satisfy following 
recurrence relation: 

00 

(e-a-u-2j)e v = £ [c^ +1 + (« + u)C% + C^ 1 ] e 2k , (3.12) 

k=j+l 



The Relativistic Linear Singular Oscillator 



8 



where C™ are binomial coefficients. From (|3.12|) it follows that the power series (|3.11|) 
will be terminated at the term e2 n {ip) 2n if condition e = e n = 2n+a + v, n — 0, 1, 2, 3, . . . 
holds. This gives following quantization rule for the energy spectrum for relativistic 
singular oscillator (j3.1|) : 

E n = hue n = hu(2n + a + u), n = 0, 1, 2, 3, . . . . (3.13) 

Hence in the case of ()3.13|) solutions of the equation ()3.1()j) coincide with the 
continuous dual Hahn polynomials 

n(p)=n n (p) = S n (p 2 ;a,v,^\ , (3.14) 



defined with relation [2Zj 

S n (x 2 ; a, b, c) = (a + b) n (a + c) n 3 F 2 



-n, a + ix, a — ix 
a + b, a + c 



(3.15) 



where (a) n = a(a + 1) • • • (a + n — 1) = T(a + n)/T(a) is Pochhammer symbol. 

Continuous dual Hahn polynomials 1)3.15)1 satisfy the three-term recurrence relation 

(A n + C n -a 2 - x 2 ) S n (x 2 ) = Aj n+1 (x 2 ) + C n S n ^ (x 2 ) , (3.16) 
where A n = {n + a + b) {n + a + c) , C n = n{n + b + c — 1) and 

* n [X > - (a + b) n (a + c) n ■ 

Hence normalized wavefunctions for the stationary states of the relativistic singular 
linear oscillator have the following form: 

MP) = c„(-p) (Q) ^ p r (u + ip) S n (p 2 ; a, u, ij , (3.17) 



r (n + a + v) T (n + a + |) T (n + v + \) n\ ' 
Wavefunctions ()3.17)) are orthonormalized as follows: 

^n(pW m {p)dp = 5 nm . 



One can verify that the difference Hamiltonian of equation ()3.3j) 

H = mc 2 [a + a~ + oo (a + v)] (3.18) 
may be factorized in terms of the operators, having the form 




Figure 1. The behaviour of the ground state wave function of the relativistic linear 
singular oscillator Ij3.17|l for values of speed of light c = oo (the non-relativistic case), 
4 and 0.25 (m = u) = % = 1). Real parts are shown in the left plots, imaginary parts 
in the right plots. 



1 
1 

71 



-ir) 



a 



e z Up - cj e 2 " (u + ip) (H 



a 



e 2 9p - u (v - ip) (1 )e2 



(3.19a) 
(3.196) 
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They are the pair of hermitian conjugate operators. Using a and a + , one can 
construct a dynamical symmetry algebra of the relativistic linear oscillator ([3.3)1 . 



non-relativistic 



c=4 



-10 





c=2 



10 
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-10 



-10 





c=l 



10 




10 



c=0.25 



10 



Figure 2. The behaviour of the ground state energy- level of the relativistic linear 
singular oscillator l|3.13fl for values of speed of light c = oo (the non-relativistic case), 4, 
2, 1, 0.5 and 0.25 (to = w = K = 1). Real parts are shown by solid lines, imaginary parts 
by dashed lines. We see that, decreasing of c changes appearance point of imaginary 
part ("collapse" point g < — g (l + ^pr)), too. 



We note that, hermiticity condition of Hamiltonian H imposes a restriction to the 
values of the quantity g. Indeed, from expression (|3.18)1 it follows that H is hermitian 
operator and its eigenvalues E n (j3.13|) are real only in the case, when a and v are real or 



complex-conjugate quantities. Therefore, g must satisfy the condition g > — 



h 2 u 



8m 



32m 
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The case g < — ("collapse") will be considered separately. 

Since the parameter \i = — > oo when c — > oo, by using (|A.1|) in Appendix 
A it is easy to show that in the limit case, when the velocity of the light approaches 
oo, wavefunctions ()3.17j) coincide with wavefunctions (jl.2j) of the non-relativistic linear 
singular oscillator. 

Energy spectrum (j3.13j) also has a correct non-relativistic limit, i.e. 

E n - mc 2 -> E™ nrel = hu{2n + d + 1) . 



By the use of (jA.2|) it can be shown that, in the g — > limit, the eigenfunctions 
(I3.17j) and the eigenvalues (j3.13)l transform into the " antisymmetrical" eigenfunctions 
and eigenvalues of the relativistic linear harmonic oscillator, considered in 

C osc (p) = C n [1/(1/ - l)V p/2 T( P + iv')p£ +1 (p; £ 

E r n elosc = hw (2n + 1 + v') , n = 0, 1, 2, 3, . . . . 
Here P%' (p; </?) are the Meixner-Pollaczek polynomials and 



' " l/ 27rAr(2n + 1 + z/') ' 



'n 




In Figured] we show the behaviour of the ground state wave function i/}q(x) (|3.17j) 
in dependence of the x and g for various values of the speed of light c. We see that 
except for the non-relativistic case, where ipo(x) is real for g > — in the relativistic 
case ground state wavefunction is complex and has a lot of number of zeros. Similarly, 
in Figure El we show the behaviour of the ground state energy- level (J3.13P in dependence 
of the g for various values of the speed of light c. From these plots we see that decreasing 
of c changes appearance point of imaginary part ("collapse" point g < — ^ — 32^7). 



4. Conclusion 



The application of the finite-difference relativistic quantum mechanics to a large class of 
physical problems requires relativistic generalizations of the exactly solvable problems 
of non-relativistic quantum mechanics. 

In this paper we construct a relativistic model of the linear singular oscillator and 
explicitly solve the corresponding finite-difference equation. We determine eigenvalues 
and eigenfunctions of the problem. They have correct non-relativistic limits. As in the 
non-relativistic case [2], the simplicity of the obtained energy spectrum suggests that a 
solution by group-theoretical methods should also be possible. Furthermore, we study 
some properties of the continuous dual Hahn polynomials. 

We hope that, the model of the relativistic linear oscillator proposed in this paper 
will be applied in future in various fields of quantum physics likewise the non-relativistic 
singular oscillator. 
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Here we present some formulas for the continuous dual Halm polynomials (|H.15j) used 
in the text. 



1. From the recurrence relation for the continuous dual Hahn polynomials (|3. 16)1 . 
it can be shown that the following limit to the Laguerre polynomials holds: 

lim ——3 n (zfji-,a,b,- \ = Ln 5 (» > (A- 1 ) 



/u->oo n!/i n \ 2 
where ao = lim a and lim (b — fi) = const. 

/i— >oo /j— »oo 



2. It can be shown that, in some particular cases continuous dual Hahn polynomials 
coincide with the Meixner-Pollaczek polynomials, i.e. 

PL + i I) = (-l) n 7^^T ^n fa; 2 ; 1, 6, ^ (A.2) 

and 



(2n+l)\ " V 2 



(A.3) 



To prove formulas (|A.2J) and (|A.3|) let's compare equations [2Z| 
[e~^ (6 + ix) e~ idx - e iip (b - ix) e idx ] y x (x) = 2i [x cos <p-(k + b) sin <p] y x (x) (A.4) 
and 

[(a + ix) (b + ix) e~ i9x - (a - ix) (b - ix) e idx ] y 2 (x) = 2ix (2n + a + b) y 2 (x) (A. 5) 

where yi(x) = P%(x,<p) are the Meixner-Pollaczek polynomials and y 2 (x) = 
S n (x 2 ; a, b, |) are the continuous dual Hahn polynomials. It is easy to verify that, 
when v 9 = f an d k — 2n + 1, Eq. (jA.4j) coincides with Eq. (jA.5|) for a = 1. It means 

An+i (*; I) = ^n^n (x 2 ; 1, 6, . (A.6) 

Comparing coefficients, for example, for x 2n+1 in the left-hand and right-hand sides 
of (TO|) we find that, 



N n = (-1Y- 



(2n+ 1)! 

One can prove relation (jA.3|) in same way. 
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3. In the non- relativistic limit, when /i = ^ > oo we have 

lim a = d H — , 
lim (V - /x) = - , 



lim (-p) (Q) = e K d +i) ln ^ , (A.7) 

lim M(p) = s/ZKe^v-^ , 



lim C„ = 1 _ p „-(„+n+i)ln» ^ 

m^oo TTA/nir (n + d + 1) 
where £ = To obtain these formulas we used the representation 

/2tt 

r(z) ~ J — e 2ln2 - 2 , 
|z| — > oo , 
for the gamma function. 



4. In the non-relativistic limit we have following limit relations for operators (j3.19aj) 
and lEHSfy 



if d + \ 
hm ^Jla =c = + £ 



lim ^fjia + = c + = -^= ( -9 e + £ - 

Hamiltonian 1)3.18)1 in this limit coincides with Hamiltonian of the non-relativistic 
singular oscillator ([1.1)1 : 

h f d 2 — - \ 

lim H = H N = hu (c + c~ + d+ 1) = — f -a 2 + £ 2 + — — -A ) . 

M ^oo 2 v i ) 
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